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Abstract 



We attempt to understand the low-mass dielectron enhancement observed by PHENIX Collab- 

cn : 

oration at Relativistic Heavy Ion Collider (RHIC) by transport peak in the spectral function. On 
the basis of the second-order formalism of relativistic dissipative hydrodynamics, we parameterize 
the spectral function in low- frequency and long- wavelength region by two transport coefficients, 
electric diffusion coefficient D and relaxation time tj, and compared our theoretical dielectron spec- 
tra with the experimental data. We study spectrum of dielectrons produced in relativistic heavy 
ion collisions by using the profile of matter evolution under full (3+l)-dimensional hydrodynamics. 
We find that the experimental data require the diffusion coefficient to be D > 2/T, with T being 
temperature. Our analysis shows that dielectrons emitted through transport process mainly come 
from high-temperature QGP phase. 

PACS numbers: 



I. INTRODUCTION 

Study of the quark-gluon plasma (QGP) is currently undertaken by Relativistic Heavy 
Ion Collider (RHIC) at BNL and by Large Hadron Collider (LHC) at CERN. At RHIC, 
signatures of the strongly correlated QGP have been accumulated through the analyses 



u 



of the collective evolution of hot QCD matter [1 - 3 ] and of the energy loss of jets and 
heavy-quarks [4], |5|]. Electromagnetic probes, i.e. photons and dileptons, are also important 
probes of the QCD matter since they are not contaminated by strong interaction and reflect 
transport property of the matter. 

Recently, PHENIX Collaboration at RHIC reported enhancement of dielectrons in low- 
mass (0.1 < m cc < 0.75 GeV) region in Au+Au collisions p, |7J. However, theoretical 
models, which were successful in reproducing the low-mass dilepton spectra at SPS (Su- 
per Proton Synchrotron), cannot explain the PHENIX data. This indicates the existence 
of yet unknown sources beyond the standard thermal radiations J8l-ll3|. By modeling the 
spectral function based on the two scenarios discussed in the literature, the dropping mass 
and the width broadening, and combining it with the full (3+l)-dimensional hydrodynamic 
evolution, we also find that our model spectral function cannot explain the low-mass en- 
hancement at PHENIX. (See Appendix lAl for details.) As pointed out in [lOj, one of the 
possible candidates which have tendency to fill the gap could be the processes considered 
in the L andau-Pomer anchuk- Migdal resummation, such as the off-shell annihilation process 
q + q + q — > q + 'y* — > g + e + + e~ [14| . However the use of perturbative picture near transition 
temperature is not necessarily justified. 

The main purpose of this paper is to study transport peak in the spectral function and 
its consequence on dilepton spectra in low mass region. Since the transport peak gives 
divergent dielectron rate in low-frequency limit u — > 0, this could be a possible source of the 
low-mass enhancement. The transport peak reflects transport property of the QCD matter 
and has not been fully considered in the above calculations. Perturbative calculations of the 
transport coefficients and spectral function in low-frequency region 15M20| require higher- 
order resummation and thus the transport peak originates from various multiple scattering 
processes. Instead of taking the perturbative result, we parametrize the transport peak 
at low-frequency and long-wavelength by a set of transport coefficients, the electric charge 



diffusion coefficient D and relaxation time tj 2l|, |22j . Then we try to constrain their values 



by analyzing the dielectron data at PHENIX with the use of the full (3+l)-dimensional 
hydrodynamics simulation and the state-of-the-art lattice equation of state. Furthermore, 
we compare the resultant constraint with the perturbative QCD estimate at weak-coupling 
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- |l8| and the AdS/CFT estimate at strong coupling 23|. We find that the main source 
of the low-mass dielectrons emitted through the transport process is the high-temperature 
QGP, not the low-temperature hadronic phase. 

In Sec. [TTl we review the basics of the relativistic hydrodynamic model. In Sec. IHfl we 
utilize the second-order formalism of relativistic dissipative hydrodynamics in the presence of 
external electromagnetic field to obtain the spectral function parametrized with two trans- 
port coefficients, D and rj. Then, we calculate dielectron production using the spectral 
function with a transport peak and compare the results with the experimental data. In 
Sec. HVt we give conclusion and outlook. 

II. RELATIVISTIC HYDRODYNAMICS 

The relativistic hydrodynamic model has been quite successful in describing collective 
flow phenomena in heavy-ion collisions at RHIC [iMsl]. Its basic equation for perfect fluids 
reads 

d^ v = 0, T^ = (e + P)uV - Pg^, (1) 

where T^ y is the energy-momentum tensor, e the energy density, P the pressure, and u^ the 
fluid velocity. The baryon chemical potential is neglected, since it is small near mid-rapidity 
at RHIC and LHC energies. 

A. Lattice equation of state and spacetime evolution 

The energy density and pressure are related through the equation of state (EoS) P = P(e). 
In the present paper, we use one of the latest EoS obtained from (2+l)-flavor lattice QCD 
simulations with a Symanzik improved gauge action and a stout-link improved staggered 
fermion action [24]. The parameterization of the trace anomaly I = e — 3P given in [24] 
reads 

I{T) { . ,. , ,, 2 , (. /o • [tanliCA • t + / 2 ) + 1] ^ (0 s 

- exp(-hi/t - h 2 /t ) • [ho + — — - 2 , (2) 
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FIG. 1: (Color online) Pressure as a function of energy 



density P = P(e) compared with the conventional 



24j. In the bag model at T = 170 MeV, pressure 



bag equation of state and the lattice equation of state 

with Nj = 3 massless free ideal QGP gas is matched to that of a hadron resonance gas including resonances 

up to A(1232). The resultant bag constant is B l ' A = 247 MeV. 

where t = T/(0.2GeV), (AoAA) = (0.1396,-0.1800,0.0350), (f , fi, f 2 ) = 
(2.76,6.79, -5.29), and (g u g 2 ) = (-0.47, 1.04). As shown in Fig. H the lattice EoS shows a 
smooth crossover from the hadronic phase to the QGP phase in contrast to the historic bag 
EoS, which has the first order phase transition. 

In solving hydrodynamic equation Eq. (CQ), we pose the initial condition at r = 0.6 fm/c 
for entropy density s and flow vector u: 



T s{r} 8 ,x±) = C ■ 9(y h - \r] s \) 

xf PP (Vs) 
u z = sinhr/s, u x 
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(3) 
(4) 



when nucleus A(B) is traveling along z-axis to negative (positive) direction. Here y^ (> 0) is 
the beam rapidity, C = 13.0 and a = 0.85 are fitting parameters, dN p ^ rt /d 2 x±,dN co \\/d 2 x± 
are defined in the Glauber model 25], and f pp (ri s ) is a parametrization of the shape of 
rapidity distribution in proton-proton (p+p) collision at ^/snn = 200 GeV. For details of 



the initial condition, see Refs. 
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29]. With this initial condition, hydrodynamic equations 



Eq. (CQ) are numerically solved assuming chemical equilibrium above freezeout temperature 
Tf, which we determine as follows. 



At T = Tf, Cooper- Frye formula 30] converts the profile of hydrodynamic evolution 



> 

CD 

C5 



10 

1 

0.1 

0.01 

0.001 

0.0001 



T f =0.15GeV(20-30%) 

PHENIX(20-30%) i — i — i 



(a) 



> 

CD 

C5 



0.5 



1 1.5 2 

p T [GeV] 



2.5 



10 



0.1 



0.01 



0.001 



0.0001 



T f =0.15GeV(40-50%) 

PHENIX(40-50%) i — i — i 



(b) 



0.5 



1 1.5 

p T [GeV] 



2.5 



FIG. 2: (Color online) Shown is the hydrodynamic calculation for the slope of p T spectra of protons with 
freezeout temperature T{ — 0.15 GeV. The p T spectra are calculated by the hydrodynamic model with lattice 
EoS. In the hydrodynamic calculation, chemical equilibrium is assumed at T > T{ = 0.15 GeV. The impact 
parameters of the collision are (a) b = 7.1 fm (20-30% centrality) and (b) b = 9.7 fm (40-50% centrality) in 
the hydrodynamic calculation. Note that overall p T spectra are scaled (by multiplying sa 1.4) to match the 
experimental data at p T = 1.0 GeV. 

(u,T) into spectra of particles and resonances (up to A(1232)) by 

d 3 p ~ (27r)3 J / 9s expb^/T f ] T l' [ } 

where di is the degeneracy factor of hadron species i, the sign in the denominator is — (+) for 
bosons (fermions), <9£ is the spacetime hypersurface satisfying T(t, rj s , x, y) = If, and da^ is 
its element. To obtain spectra of stable particles, the decay products from A(1232) are also 
included. We determine the thermal freezeout temperature 7f by fitting the experimental 
slope of the proton p T spectra as shown in Fig. [2j The slopes can be fitted well by Tf = 0.15 
GeV. The impact parameters of the simulation are b = 7.1 fm (20-30% centrality) and 
b = 9.7 fm (40-50% centrality). Note that we have multiplied a factor w 1.4 to fit the 
absolute magnitude of the experimental data at p T = 1.0 GeV. This slight difference in the 
magnitude between hydrodynamical calculation and experimental data is possibly due to an 
absence of chemical freezeout mechanism 27J in the lattice equation of state. 



B. Dielectron spectrum 

Here we briefly summarize the formula for the dielectron production from expanding 
medium in relativistic heavy-ion collisions. Since the invariant mass of our interest is 



\f(f > 0.1 GeV, we neglect the electron mass in the following. The (3+l)-dimensional 
hydrodynamic model gives each spacetime point a local flow vector u^ L (x) and a local tem- 



perature T(x). 



with the rate 
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ach spacetime volume under hydrodynamic evolution emits lepton pairs 

m 



d 3 Pl d 3 p 2 d 4 x " 2vrV exp(g°/T)-l R W ' h l ' 

Gg(q;T) = jd 4 xe^9(x°)([J^x),r(0)]) T , (7) 

J"{x) = ^u-~d 7 ^d-^s, (8) 

when observed in the local fluid rest frame. In the laboratory frame, momentum distribution 
of the lepton pairs emitted from the thermal medium in heavy ion collisions is expressed by 
the following integration: 

E l E 2 dN f . a 2 V\vl+& x -^9 ilv 



7 E 2A~4 exp(gO/T(x)) - 1 imG aU,JWJ. [J) 



d 3 pid 3 p 2 

E indicates the spacetime region of the hot QCD medium with T > Tf which is considered 
to radiate dielectrons. The A^ v " denotes a Lorentz tensor in the local rest frame boosted 
by u^ from that in the laboratory frame A 1 "*'". 

In the PHENIX experiment, with which we will compare our results later, the acceptance 
of each electron or positron depends on the detector geometry and on the particle kinematics 
pi, LZj] - We introduce an acceptance function A(p 1 ,p 2 ) to pick up the phase space where both 
the electron and positron are in the acceptance. The dielectron invariant mass spectrum 
from the expanding medium is therefore given by 

dN „ f ., , q dN 



dm ee 



2m ee Jd 3 Pl d 3 p 2 j 3 ^ 3 A(p u p 2 )6(q 2 - m 2 ee ). (10) 



III. MODEL OF LOW-MASS DILEPTONS: TRANSPORT PEAK 

One of the promising non-perturbative approaches to explain the low-mass enhancement 
is to utilize the transport theory, which parameterizes the spectral function in low-energy 
and long-wavelength region with transport coefficients. 



A. Linearized hydrodynamics with external field 

Here we construct retarded correlator G^(q;T) and corresponding spectral function 
(transport-SPF) from the linear analysis of transport equations. According to the linear re- 
sponse theory, the retarded correlator relates small perturbation 5H(t) = f d 3 xJ^(x)SA^(x) 
and system response {8J^(q))T in the linear order [211 ]: 

(8J^q)) T = -G^(q;T)8A v (q), (11) 

where f(q) = J d 4 xe tqx f(x). Hereafter we calculate the linear relation between 5J^ and 8A V 
on the basis of relativistic dissipative hydrodynamics. 

In the presence of external field 5 A 11 , conservation laws for the energy- momentum tensor 
T^ v and the electric current J M are modified to 

d u T^ = F^J V , (12) 

<V M = 0, (13) 

where we define F^ u = d^SAj, — d v 8A^ E = — V5A° — ^, and BeVx 8 A. In relativistic 



viscous hydrodynamics in Landau frame 33j, |34|, the energy- momentum tensor T^ u and the 
electric current J M are decomposed as 

T"" = eu^u u -{P + U)A^ + n^, (14) 

jv = pu » + i/\ (15) 

A"" = sT - Wu u , (16) 

with bulk pressure n, shear stress tensor tt ijj/ , and dissipative electric current z/ M satisfying 
ix^Uv = 0, ti£ = 0, v^u^ = 0. The scalar quantities e, P, and p satisfy the equation of state 
P = P(e,p) derived in equilibrium. The entropy current in the second-order formalism 



constructed by Israel and Stewart [35|, |36j is decomposed as 



s» = su» - ^V - -(a Ui^ + a x ^ v v v ) - —%Ii 2 - fti/V„ + /WO, (17) 



with positive transport coefficients a a (a = 0, 1) and /?& {b = 0, 1, 2) and chemical potential 
p for the electric charge. Here s is the entropy density s(e, p) derived from the equation of 
state. Calculating the divergence of the entropy current up to quadratic order in n, n^ , z/ M , 



and F^ u and in derivatives of e, p, and -u M , we obtain 
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(18) 



where / = u^d^f . The second law of thermodynamics requires the following constitutive 
equations: 



-n = c(dX + ao0X + A>n), 

7i> = 27]{{d ll u v - aid^ -/3 2 7J>)) 



i/* = ff A w 



Td p ( ^J + F pct m ct + a <9 p IT + ai 9 ff ^ - Piv P 



(19) 
(20) 
(21) 



with bulk and shear viscosities (, i] (> 0) and electrical conductivity a (> 0). Here ((B^)) 
stands for a spatial, symmetric, and traceless tensor extracted from a general tensor B^ v : 

B pa + B ap A p(7 A al3 B a p 



((B^)) =A w A l 



(22) 



When the system is close to equilibrium, the energy density e(x), the charge density p(x) 
and the flow vector u^(x) slightly deviate from their equilibrium quantities: 



e(x) = e + Se(x), 
p(x) = p + Sp(x), 



uHx) 



1, 5u{x)). 



(23) 
(24) 
(25) 



Neglecting for simplicity the couplings «o,i between different dissipative modes, the dissipa- 
tive part of the electric current v^{x) is then given by 

(26) 

(27) 



u°(x) = 0, 






u(x) = —a 


[TV(|) 


-E + ^dtU 



in the linear order in 5e, Sp, 5u, and A**. Explicitly in terms of 8e, and Sp, 

-//i\ 1 (de ude\ - 1 f dp p dp\ -. 

TV (t) = x{df + Td~p) V6p -X{df + Td-p) We ' 



X = 



de dp de dp 
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(28) 
(29) 



In the quark-gluon plasma with vanishing baryon density, p = |^ = -^ = holds and the 
electric current is obtained as 

J»( x ) = (8p(x),v(x)), (30) 

-> -> dv 

V{x) = aE-DV5p-rj — , (31) 

D = -, rj = fra, X = ^~- (32) 

X dp 

The dynamics of the electric current J^(x) is obtained in a closed form, while with finite 

electric charge density there arises a coupling between the sound mode 8e(x) and the dis- 

sipative electric current v^(x) even after ignoring 0:0,1- This is because energy current, or 

momentum, carries electric charge at finite density p 7^ 0. 

Combining Eqs. f fTBl . f|T5|) . and f l31|) . we obtain 

jo,, = xD{-k 2 5A {q)-ujk^A l {q)) 
—Tju 2 — iio + Dk 2 



-tjcg 2 — zcg + DA; 2 V ^ 2 



+— : -MM + yf^ > ( 34 ) 

-tjcg - z \ k 2 ) 

where g^ = (w, fc) and upper and lower Lorentz indices matter, k x = k t=l = —ki=\ for 
example. The retarded correlators are given by 

rrk 2 

Gg( W , S; T) = G»(., S; T) = _^ 2 ^ +gtf , (36) 

Gg( U , £; T) = -J=^ (* - *£) + ^ ** (37) 

— TjCG -! \ AT / — TjW Z — 2CG + iJAr At 

and we obtain two independent components of the spectral function 

MG<%; T) s -gtaGgfe T) = --^-^L—, (38 ) 

ImG<, T) («; T) = \ (lmG£>; T) - ImG^'fa; T)) = -j^y- (39) 

Here we have demonstrated our calculation in one-flavor case, but generalization to multi- 
flavor case can be achieved straightforwardly. We discuss in Appendix [B] that the transport 
equation Eq. (13T1) holds in 3-flavor case under three assumptions: (i) There is no coupling 
between the dissipative modes in the decomposition of the entropy current, (ii) Transport 
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FIG. 3: Scaled transport-SPF <j(oj, 0; T). It diverges as a oc 1/uj in the limit u — >• 0. 

coefficients are flavor-independent, and (iii) Quark number susceptibility matrix is propor- 
tional to the unit matrix Xij oc 5^. Here Xij — f^ an d n\ (/Xj) is a quark number density 
(quark number chemical potential) for a flavor i. Xij — f^ 1 1S measured by a lattice QCD 
simulation [37| and is found to satisfy Xij ^ $ij approximately 

Let us consider the scaled spectral function cr(q;T) = — haG^ l (q; T)/3q 2 for the 



transport-SPF. At k = 0, 



a(u,0;T) 



XD 



2cu tW + 1 : 



(40) 



and the scaled spectral function becomes larger for larger \ an d D and for smaller tj. This 
can be physically understood as follows: The number of photons and dileptons is large when 
there is strong electric current. When D is large and tj is small, strong electric current is 
swiftly induced by the electric charge fluctuation, which is characterized by the susceptibility 
X- As shown in Fig. [31 the scaled spectral function shows divergence a oc u~ x at u ~ 0. 

B. Diffusion coefficient and relaxation time 



Our spectral function is parameterized by three unknown functions of temperature 
x(T), D(T), and tj(T). The electric charge susceptibility x can be expressed in terms 
of the quark number susceptibility by 

2 



_ 4 1 1 

X ~ gXuU + gXdd + ~Xb 



4 4 2 

gXud - -Xus + -Xds 



:Xv 



which can be fitted to the result of lattice QCD simulation [37| as follows 



X (T,/i = 0)^0.28T 2 



1 + tanh 



0.15T* 



(41) 
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We take T* to be the temperature at the center of the chiral crossover transition T* 
GeV (see Fig. 4 (right) in 38|). We parameterize D and tj as 



D OC -,Tj oc -, 



0.155 



(43) 



for dimensional reason [39]. To see the typical magnitude of the transport coefficients of 
weakly and strongly coupled plasmas, we refer to the results obtained by perturbative QCD 
id-llal with g « 2 

0.150 
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(44) 
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and those by AdS/CFT [23 1 
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Here D^ CD stands for quark number diffusion coefficient. We can show D = D q within the 
same approximation as used in the derivation of the constitutive equation in multi-flavor 
case. (See Appendix [Bl) 



C. Comparison to experimental data 



The minimum bias theoretical spectrum is calculated by 
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(47) 



neglecting peripheral collisions with centrality 60-100 % because their contributions are 
small. 

In Fig. H] (a)-(d), we compare the dielectron spectrum dN/dm ee calculated with the 
transport-SPF with the experimental data. The parameterization for the diffusion coef- 
ficient is (a) D = 10/T, (b) D = 5/T, (c) D = 2/T, and (d) D = 1/T and that for the 
relaxation time is varied in 7j = (0.0-1.0)/T. We add the dielectrons from hadronic decays 
after freezeout (hadronic cocktail) [7] to the thermal contribution from the transport-SPF. 
We find that our results with (DT,tjT) = (10,0.5), (5,0.2), and (2,0.0-0.1) give dielectron 



11 




o 

E 



0.2 0.4 0.6 0.8 
m ee [GeV] 




0.2 0.4 0.6 0.8 
m ee [GeV] 



DT=2, TjT=0.0 

TjT=0.1 

Cocktail 

PHENIX('IO) 




DT=1,TjT=0.0 

TjT=0.1 

Cocktail 

PHENIX('IO) 



0.2 0.4 0.6 O.i 
m ee [GeV] 




0.2 0.4 0.6 0.8 
m ee [GeV] 



FIG. 4: (Color online) We compare in (a)-(d) the dielectron spectrum dN/dm cc calculated with the 
transport-SPF and the corresponding experimental data jj. In the transport-SPF, diffusion coefficient is 
fixed to (a) D = 10/T, (b) 5/T, (c) 2/T, and (d) 1/T. The relaxation time tj is varied in (0.0-1.0)/T. We 
take into account the contributions from the hadronic decays after freezeout (denoted as "Cocktail" in these 
figures). We only plot the statistical errors for experimental data. 

spectra comparable to the experimental data. These values are rather different from the 
theoretical estimates, (L>PQ CD r, Tj pQCD T) « (4, 15) and (D AdS / CTT T, Tj AdS / CFT T ) = (i |2). 
From (c) and (d), we also find that experimental data set lower bound for the diffusion 
coefficient D > 2/T because the dielectron rate is largest with tj = for fixed D. 

Next we study the dielectron spectrum in more detail. In Fig. |5] (a) and (b), we show 
the dielectron invariant mass spectrum with the dielectron transverse momentum restricted 
in (a) < p T < 0.5 GeV and (b) 0.5 < p T < 1.0 GeV. We calculate these spectra with 
parameter sets (DT,tjT) = (10,0.5), (5,0.2), and (2,0.0-0.1) in the transport-SPF. Clearly 
all of our calculations for low p T -window (a) undershoot the experimental data in 0.3 < 
m ee < 0.7 GeV, while for high p T -window (b) they all show good agreement with the data. 

In Fig. El we show the dielectron transverse mass (m T ) distribution at mid-rapidity 
y p = 0. The calculation is performed with (DT, rjT) = (10,0.5), (5,0.2), and (2,0.0-0.1). 
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FIG. 5: (Color online) In (a) and (b), we compare our calculation of the dielectron spectrum with restricted 
p T -windows and the experimental data 7j. The dielectron spectrum dN/dm cc is restricted in (a) 0.0 < p T < 
0.5 GeV and (b) 0.5 < p T < 1.0 GeV. We calculate with parameter sets (DT,t 3 T) = (10,0.5), (5,0.2), and 
(2,0.0-0.1) in the transport-SPF. The hadronic decays after freezeout are taken into account (denoted as 
" Cocktail" in these figures) . In (a) and (b) , we only plot the statistical errors for experimental data. 
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FIG. 6: (Color online) We compare the transverse mass spectrum calculated with the transport-SPF and 
the corresponding experimental data [7]. The parameter sets (DT,tjT) are the same with those in Fig. [5] 
In the experimental data, other sources such as hadron decays are subtracted and PHENIX acceptance is 
corrected. We only plot the statistical errors for experimental data. 

We first calculate the dielectron p T -spectrum by integrating in the invariant mass range 
0.3 < m cc < 0.75 GeV for each centrality. We then divide the dielectron p T -spectra by half 
of the number of participant nucleons (N part /2) for each centrality. We take the horizontal 
axis to be m T — m = \ p\ + vrc^ — mo to obtain m T -spectra. Here mo is the mean value 
of the dielectron invariant mass with minimum bias mo = (m ee ) in 0.3 < m ee < 0.75 GeV, 
which turns out to be mo ~ 0.47 GeV in our calculation. Experimental data are also 
obtained by the same procedure after subtracting other sources such as hadron decays and 
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FIG. 7: (Color online) Plotted is dielectron production from the matter with T > 0.3 GeV, 0.2 < T < 0.3 
GeV, and T < 0.2 GeV. The parameter for transport-SPF is (DT, t 3 T) = (10, 0.5). The PHENIX acceptance 
is taken into account. 

correcting PHENIX acceptance. The experimental data exceed the theoretical spectra at 
m T —m Q < 0.3 GeV while the former are comparable with the latter at m T — m > 0.3 GeV. 
This is consistent with what we found in Fig. [5j 

In Fig. [7J we show dielectron production from the hot matter with different temperature 
range: T > 0.3 GeV, 0.2 < T < 0.3 GeV, and T < 0.2 GeV. Clearly the intermediate- 
temperature QGP phase with 0.2 < T < 0.3 GeV is the main source for the dielectron 
production in spite of its small spacetime volume. This makes a sharp contrast with the 
expectation that the low-mass dielectrons are radiated from the low-temperature hadronic 
phase due to the larger spacetime volume [8|] . Furthermore, this tendency may be consistent 
with the previous analyses at SPS, where the transport peak does not play an important 
role. 



IV. CONCLUSION AND OUTLOOK 

In this paper, we have studied the dilepton production from the hot dynamical medium 
in relativistic heavy ion collisions. In particular, we have explored the novel relation between 
the low-mass dielectron enhancement and the transport property of the hot QCD matter. 
Because the dielectron rate from the transport-SPF is divergent at u ~ 0, it is reasonable 
to expect that the transport property is relevant to the low- mass dielectron enhancement. 
Also the transport theoretical approach enables us to treat the problem non-perturbatively. 

We have conducted our analysis in the following way. We start from the second-order 
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FIG. 8: (Color online) We plot parameter sets (DT, rjT) with which the experimental data dM/dm ee 
are consistent. We also plot two limiting sets obtained by weak-coupling (pQCD) and strong-coupling 
(AdS/CFT) calculation. 

formalism of relativistic dissipative hydrodynamics, which introduces diffusion coefficient D 
and relaxation time tj. By performing linear analysis, we parameterized the low-frequency 
and long- wavelength region of the spectral function with these transport coefficients. We cal- 
culated the dielectron yield by combining the transport-SPF and the full (3+l)-dimensional 
hydrodynamic medium evolution with the lattice EoS and attempted to extract the trans- 
port coefficients from the data. We summarize our results in Fig. [HJ We found that trans- 
port peak with parameter sets (DT,tjT) = (10,0.5), (5,0.2), (2,0.0-0.1) can reproduce 
the experimental data dN/dm ee . These obtained parameter sets are rather different from 
both weak-coupling pQCD calculation and strong-coupling calculation on the basis of the 
AdS/CFT correspondence. We also showed that in order to explain the data, the diffusion 
coefficient must be D > 2/T. While these parameter sets could reproduce the more detailed 
experimental data dN/dm ee restricted in 0.5 < p T < 1.0 GeV and (l/27rm T )dN/dm T dy p 
in m T — m > 0.3 GeV, they could not reproduce those restricted in 0.0 < p T < 0.5 GeV 
and in m T — m® < 0.3 GeV. In view of these calculations, we conclude that the low-mass 
dielectrons at RHIC has not been fully understood theoretically despite the fact that the 
transport peak of the spectral function has a right tendency to enhance low-mass dileptons. 
Meanwhile we also found that a large portion of the thermal dielectron radiation comes 
from the high-temperature QGP phase with T > 0.2 GeV because the fluctuation of electric 
charge is larger at higher temperature. Note that RHIC has explored the temperature 
range T > 0.2 GeV for the first time ever. This tendency of dielectron production from 
the transport peak may be consistent with the dielectron production at SPS, which can be 
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explained without the transport peak. 

We list several future implications of interest: (i) First-principle calculation of transport 
coefficients and spectral function by lattice QCD is highly desirable. Calculation of the 
spectral function with finite momentum may give a clue to understand the origin of low- 



mass die 
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ectron enhancement. Some recent developments in such a direction can be seen in 



50|. (ii) The event-by-event fluctuation of the initial geometry, which is expected 



to remain large at the early QGP phase, would enhance the dielectron production from the 



transport peak, as it does for real photon emission from high temperature region [5l|]. (iii) 
Since high-temperature QGP phase emits a large portion of thermal dielectrons in spite of 
its small spacetime volume, the dielectron production from pre-equilibrium stage, or the 
so-called glasma, might also be important and thus needs to be evaluated. 
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Appendix A: Model of low-mass dileptons: vector meson dominance 

Here we assume vector meson dominance on the hadronic coupling to virtual photon 
and study its consequence on the dilepton spectra. Two standard scenarios, the dropping 
mass and the width broadening of the vector meson spectral function (VMD-SPF) will be 
considered. 
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1. Vector meson at finite temperature 



We define the dimensionless spectral function at finite T as a(q; T) = —ImG^(q; T)/3q 2 . 
In the vacuum at T = 0, it can be is measured by the cross section of the electron-positron 
annihilation (e + e~ — > hadrons) and can be decomposed as 

a{q- T = 0) = a p {q- T = 0) + i<r w (g; T = 0) + ^(g; T = 0), (Al) 

where a p (q; T), cx w (g; T), cr<^(g; T) and <r(g; T) correspond to the current correlation of J^ = 
i(u7^u - d7^d), J£ = |(u7^u + dyd), J£ = s^s and J^ = J» + \j£ - §J£, respectively. 
The spectral function at zero temperature, a(g; T = 0), has hadronic contribution at low 
q and perturbative QCD continuum at high q with an approximate continuum threshold at 
1 — 1.5 GeV and may be parametrized as [4( 



* v (q;T = 0) 



fv r v m v 



v <q < 8 1V , 



(A2) 



IV 



< q" 



with V = p,u,(p. Here s ov corresponds to the threshold associated with the decay into 
two pions (280 MeV) or three pions (420 MeV), while s iv corresponds to the continuum 
threshold. We take the height of the continuum cy from the leading order perturbation 
theory and fit the other parameters (the resonance mass m v , the resonance width T v , the 
resonance height f v and the continuum threshold s iv ) by using experimental data in 41 ]. 
These parameters are listed in the Table HI 

Since the p channel is a dominant source of the thermal component of the low-mass 
dileptons and is expected to receive the medium effect most strongly, we hereafter restrict our 
analysis on the medium modifications of <J p (q; T) at finite T within the same parametrization 
as Eq. flA2j) . Then the simplified form of the dilepton production rate reads 



EiE2uR i 



i+i- 
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d 3 pid 3 p2 



6ttV 
2 r 
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/ BE (g°;T) 
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2vr 4 



1 + 



2m\ 
q 2 



<t p (q;T) 



/be(<z°;T), 



(A3) 



where we implicitly assume that the transverse and longitudinal spectral functions are the 
same even at finite T. 

Let us now consider so-called dropping mass scenario where the mass parameters, m p (T) 
and s 1 (T), are assumed to scale with the chiral condensate 42M45I] : 



m P (T) 
m p (0) 



s i P ( T ) __ (qq), 



s lp (0) 
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TABLE I: The parameters of the vacuum spectral functions given in the unit of GeV except for dimen- 
sionless parameters c v . They are obtained by fitting the experimental data 41]. 
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FIG. 9: (Color online) Shown is the scaled spectral function <j p (uj, 0; T) as a function of u at several different 

temperatures. The parameters in the spectral function are (T*,A,T p ) ~ (0.155 GeV, 0.025 GeV, 0.3 GeV). 

The scaling function in Eq. (IA4I) is obtained by fitting the latest lattice data (Fig. 4 (right) 
in [381) by the following ansatz with (T*, A) = (0.155 GeV, 0.025 GeV): 






1 — tanh 



c 

A 



(A5) 



We assume that the pion mass does not change appreciably below T*, and adopt the 
following prescription for the low-energy threshold, s 0p = min.{2rra 7r , m p (T)}. Instead of 
introducing a parametrization for the T-dependence of T p (T), we vary its value in the range 
0.15-0.45 GeV to see its effect on the dilepton yield. Once it is given, remaining parameter, 
f p (T), is constrained by the QCD spectral sum rule [461 ]: 



/■oo 

/ dtu 2 [a p (tu,6;T)-c p ] = 0. 
Jo 



(A6) 



Shown in Fig. [9] is a resultant spectral function in the p channel with a parameter set, 
(T*, A,T P ) = (0.155 GeV, 0.025 GeV, 0.3 GeV). 
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FIG. 10: (Color online) Shown in (a) and (b) are sum of the dielectron spectra from the thermal medium 
and the contributions from the hadronic decays after freezeout. The latter is denoted as "Cocktail" in these 
figures. The parameters of the VMD-SPF arc T p = 0.15,0.3,0.45 GeV and (a) m p (T) oc (qq) T (dropping 
mass) and (b) m p (T) = 0.77 GeV (collisional broadening). These spectra are compared with experimental 
data of dielectron spectrum with minimum bias |7j. In (a) and (b), we only plot the statistical errors for 
experimental data. 



2. Dielectron spectra 

In Fig. [TQl we compare the theoretical spectrum with the experimental data by taking 
into account the contributions from hadronic decays after freezeout (hadronic cocktail) given 
in [7|. The PHENIX acceptance is already taken into account. We show in Fig. [10] (a) the 
thermal dielectron spectrum for T p = 0.15, 0.3, 0.45 GeV and m p (T) = 0.77 GeV (collisional 
broadening). It is clear from this figure that the dielectron emission with collisional broad- 
ening only cannot explain the experimental data in the low invariant mass region m cc < 0.6 
GeV. We show in Fig. [TD] (b) the thermal dielectron spectrum for T p = 0.15,0.3,0.45 GeV 
and m p {T) oc {qq)x- Even with the dropping mass, the dilepton yield undershoots the ex- 
perimental data substantially in the low invariant mass region m ee < 0.6 GeV. The above 
findings are consistent with the previous attempts to reproduce the low-mass dileptons at 



PHENIX 
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Appendix B: Extension of Transport Spectral Function to Multi-Flavor Case 

We derive the transport spectral function for multi-flavor case. In this derivation, it is 
more convenient to treat flavor current for each quark species N^. We start from conser- 
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vation laws in external field 5 A 11 which couples to each flavor current and gives perturbed 
Hamiltonian 6H(t) = J d 3 xJ2f qfN tM (x)SA IJ ,(x): 



f 
d^ = (/ = u,d,s) 



(Bl) 
(B2) 



with Ffj, u = d^5A u — d u 5A^ and a» = 2/3, gd,s = —1/3. According to relativistic viscous 
hydrodynamics in Landau frame 33], we decompose T^ v and N't as 



T"" = euV -(PH-n)A^ + 7i 



\1V 



n? 



A 



IM' 






(B3) 
(B4) 
(B5) 



with bulk pressure II, shear stress tensor tt' j ' u , and dissipative flavor current uf satisfying 
-n iiV u v = 0, 7r^ = 0, u^u^ = 0. The entropy current in the second-order formalism 23l. l35l. |36| 
is decomposed as 



s M = su» 



^ 



/ 



u 
2T 



f ( ^n 2 - E A r tyr# + /^% J - \ ( E «oM + E "I*'""!,* )(B6) 

with coupling coefficients a Q ,a[ and /3 0i 2 (> 0), a positive semi-definite matrix /3[ , and a 
chemical potential [i$ for each flavor. Divergence of the entropy current by using Eqs. (IB3P 
and flB3p up to second order deviation from equilibrium is 

— f A^d^Uu + /3 n + E a o<V/ I + ~7fr ( 9^u„ - fair,* - E a i d n v ft' 



d^ = - 



-E 



1L 
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TO, (^-) + qfF^u" - E P*/*?* + «o^B + a{d u n; 

f 



(B7) 
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where / = u^d^f . Constitutive equations which ensure the second-law of thermodynamics 
are obtained as follows 



n = c ( A""^ + /3 tl + J2 <4drf j , 



ir^ = 277 ( ( &*u v - (3 2 ti^ - J2 <A^v\ 



/ / /' 



(B8) 
(B9) 



*y = '/Assy 



T0' (^) + g/^X - ^ /3 1 / ' / "/>;» + a f 'dm + a(d^ p (BIO) 



with bulk and shear viscosities (, rj (> 0), flavor conductivity <j/ (> 0), and a positive 
semi-definite flavor mixing matrix k" . ((B^)) stands for 



((S^)) = A w A l 



2 3 

We perform linear analysis in terms of Se, Srif, and Su defined as 



(Bll) 



e(x) = e + 5e(x), 
n A x ) = n f + 6rif(x), 
u»(x) = (1, Su(x)), 



(B12) 
(B13) 
(B14) 



and external vector field SA P . For simplicity we neglect the couplings between different 
dissipative terms: a^ = a{ = 0, f3( = /3^Sff, and n*? = n^Sff (assumption (i)). The 
dissipative flavor current Vf(x) is then given by 



u%x) = 0, 

V f {x) = -<Tf 



TV (^) - q f E + (5[d t u f 



(B15) 
(B16) 



in the linear order in Se, Snt, Su, and A p . We define susceptibility matrices X and Xf as 
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(B17) 
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Then the dissipative flavor current z/f is obtained as 



ft 



0, u f = q f a f E - p fjl a f d t u f - a f [~^,{e f }j X 1 



( VSe\ 

Wn d 
\V5n s J 



(B18) 



R} = {1,0,0}, {e d } = {0,1,0}, {ej = {0,0,1}. (B19) 

We make another simplification that Gf and (5[ are flavor independent (assumption (ii)): 

a f = a, P(=Pi, (B20) 

and restrict ourselves to the situation with vanishing flavor chemical potential \x$ = 0, which 
decouples flavor dissipation and sound mode propagation. We then arrive at 







(V5n\ 



Vf + P^jfff = 1S aE ~ a { g /} X t x 



(B21) 



\VSn B j 

So far we have simplified transport coefficients in order to reduce unknown parameters 
and to avoid mode couplings. In order to obtain simple constitutive equation for diffusive 
electric current in multi-flavor case, it is necessary to simplify thermodynamic quantities, 
namely susceptibility matrix to be Xf « (dn n /dfi n ) 1 = x uu l (assumption (iii)). This yields 



r\ — 

Uf + f^icr—Uf = qfaE VSrif. 



(B22) 



Note that by this simplification we also obtain electric charge susceptibility x m terms of u- 
quark number susceptibility x — (2/3)% uu . Then constitutive equation for diffusive electric 
current is given by 



Sn = y^ j qfSnf, v = y^qfisf, V+ fiia^v = -aE — -^V5n. 

f f 
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Comparing with Eq. ( 1HT|) . transport coefficients for the electric current are related to those 
for flavor currents by 



rj = p x a, a = -a, D = — , 
3 3x 



(B24) 



satisfying a = xD. 
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We are also interested in the constitutive equation for diffusive quark number current. 
Then we need to modify Eq. (IB23I) by substituting qf — > 1, 5A^ — > 5A^, and E —¥ E q , where 
5A^. denotes external vector field that couples with quark number current. Note also that 
total quark number susceptibility is Xq = 3%uu because of the simplification for Xf. Then 
we get 

r\ Q — 

5n q = ^ <W> ^q - ^2 ^' ^ + P 1& ~&^ = 3 ^ ^ 5n v ( B25 ) 

/ / at Xq 

r q = fca, a q = 3a, D q = — , (B26) 

Xq 

satisfying a q = x q .D q . Therefore we can verify D = D q , which we use in the text. 
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